TUM/T39-99-20 
RUB-TPII-11/99 



NLO corrections and contribution of a tensor 
gluon operator to the process 7*7 nn 



bA 



N. Kivel[|"'^ L. Mankiewicz"'^ M.V. PolyakoV 



""Physik Department, Technische Universitdt Miinchen, D-85747 Garching, 

Germany 

''Petersburg Nuclear Physics Institute, 188350, Gatchina, Russia 
N. Gopernicus Astronomical Genter, ul. Bartycka 18, PL-00-716 

Warsaw, Poland 

Institut fiir Theoretische Physik 11, Ruhr- Universitdt Bochum, D-44^80 

Bochum, Germany 



February 1, 2008 
Abstract 

We have calculated the NLO corrections for the leading twist amplitude of the 
reaction 7*7 tttt at large virtuality of the photon and at W"^ <^ Q^. With a 
simple model for the two-pion quark and gluon distribution amplitudes we have 
estimated that the NLO effects may reduce the leading-order amplitude by as much 
as 30%. We have found that besides the usual twist-2 quark and gluon operators 
the NLO amplitude receives contribution from a distribution associated with twist-2 
tensor gluon operator. The new gluon distribution can be investigated by studying 
the azimuthal asymmetry of the cross section for e + 7 — > e' + vrvr. 
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Introduction 

Investigation of the reaction 7*7 tttt at large photon virtuahty but 
small invariant mass of the two-pion system [W"^ <C Q^) opens a new way to 
probe partonic structure of the pion ^. Under these conditions the cor- 
responding amplitude can be factorized ^ into perturbatively calculable 
coefficient functions and nonperturbative soft matrix elements - two-pion dis- 
tribution amplitudes (27rDA). These nonperturbative matrix elements con- 
tain rich information about the partonic structure of the pion. 

An interesting feature of the reaction 7*7 tttt is that it is related 
by crossing to the deeply virtual Compton scattering (DVCS) on a pion 
target. In particular, the 27rDA's which occur in the description of the 7*7 — > 
TTTT leading-twist amplitude can be related to the corresponding quark and 
gluon skewed distributions in the pion @, 0. As we shall discuss below, 
the coefficient functions in the reactions 7*7 —>■ tttt and 7*7? — > 77? are also 
related to each other by crossing. 

From another point of view, the reaction 7*7 — * tttt can be considered 
as a minimal generalization of the transition form factor 7*7 tt^. Due 
to the different isospin of the final hadronic state these two processes are 
complementary to each other - the former probes isosinglet quark and gluon 
structure of the pion while the latter is sensitive to the quarks in the isovector 
state. 

The lowest order QCD contribution to the reaction 7*7 tttt has been 
recently analyzed in Ref. In this Letter we extend this analysis to the 
next to leading order (NLO) in asiQ"^)- The leading-order, leading-twist 
amplitude accommodates only collisions of two photons with equal helicities, 
related by crossing to the photon helicity-conserving DVCS on a pion. At 
the NLO there is a new contribution from collisions of photons with opposite 
helicities, related to photon heli city-flip contribution to DVCS. In terms of 
the Operator Product Expansion this amplitude singles out a twist-2 tensor 
gluon operator which cannot be studied in deep-inelastic scattering (DIS) 
on a nucleon target. It was introduced by Jaffe and Manohar in their 
studies of DIS on a spin-1 target. Its relation to DVCS on a nucleon and 
to the current fragmentation in DIS have been considered in Refs. |Q and 
0, respectively. As we discuss below it can be also studied via azimuthal 
asymmetries in the hard photoproduction of a pion pair. 

General definitions 

Reaction 7*7 tttt is schematically depicted in Fig.l. The kinematics 
can be conveniently described in terms of a pair of light-like vectors n, n* 
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which obey 

^2 = = 0, n-n* = 2, 

n = (l,0,0,-l),n* = (l,0,0,l) ^ ^ 

and define longitudinal directions. Vectors n, n* can be also used to define 
the metrics tensor in the transverse space 

\ n ■ n* J 

such that {-g^'')±ni, = (-^^^)±n* = 0. 

l(q) T^(ki) 



Kq) Tt(k2) 
Fig. 1 

Initial and final states momenta can be conveniently decomposed as 

-Q . 2 n2 ' (Q '2 n 

q^^ri + -n ,q = -Q ^ ^ 2Q'^' ? = 0' 

i^ = ? + ?' = |rz* + ^n, P^ = W^ (3) 

Here we have introduced the total momentum P of the pion pair, P — ki+k2, 
and the momentum fractions ( = {ki-n)/{P ■n)X = {k2-n) /{P ■n)j( + ^ — 1, 
which describe the distribution of longitudinal momentum between two pions. 
Alternatively, 

2C - 1 = /3 cos 6cm , 

where ^cm is the polar angle of the pion momentum in the CM frame with 
respect to the direction of the total momentum P and /3 is the velocity of 
produced pions in centre of mass frame 



P 



4m2 
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The amplitude of hard photoproduction of two pions is defined by the 
following matrix element between vacuum and two pions state 



T^'^ = i d*xe-"^27r{P)\TJ''{x/2)r{-x/2)\0) , q = -{q + q') 



(4) 



where J^{x) denotes the quark electromagnetic current. Hard photoproduc- 
tion corresponds to the limit ^ ^qcd; where the amplitude (|) can 
be represented as an expansion in terms of powers of 1/Q. According to 
the factorization theorem [|l], ^ the leading twist term in the expansion can 
be written as a convolution of hard and soft blocks, see Fig. 2. The coeffi- 
cient functions, defined by hard blocks, can be calculated from appropriate 
partonic subprocesses 7* + 7— i>g-|-gor7* + 7— ^^f-l-^f. 




Soil 



+ power corrections 



Fig. 2 Factorization of hard and soft contributions to the amplitude (^. 

Soft blocks represent two-pion distribution amplitudes (27rDA), or prob- 
ability amplitudes for a conversion of a pair of collinear partons (quarks or 
gluons) into the final pion pair. The quark 27rDA is defined as a matrix 
element of a light-cone quark string operator: 

(7r7r(P)|g(x)ng(0)|0),2=o = P ■ n f' dz <l>^(^, C, W^) e^^^'^ , (5) 

Jo 

with n = 7^n^. Here and below the appropriate path-ordered exponential 
between points x and is tacitely assumed. Due to the charge-conjugation 
invariance the amplitude of the process 7*7 — > vrvr depends only on the quark 
27rDA with pions in isospin zero state: 

q 

Beyond the leading order, pions can originate from a pair of collinear 
gluons as well. Here two different gluon 2'7rDA's enter. The distribution 
amplitude $'^(2;, (, W^) defined as 

(7r7r(P)|n'^n'^G^^^(x)G^"'^(0)|0)^2=o = {n ■ Pf dz <!}^{z,C,W^) e'^^^'^^ , 

(7) 
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describes the situation when pions arise from two gluons with opposite he- 
hcities, see below. As we shall discuss in the following, it can be related by 
crossing to the ordinary gluon distribution in a pion. The second distribution 
amplitude (, W^) can be defined as 

(ttttIP) Kn'' G^^^ {x)G%^ (0) |0)..=o = 

Jo 

Here the indices a±j3± are transverse with respect to vectors n, n* and occur 
in a symmetric, traceless combination. represents fragmentation of two 
gluons with equal helicities into the pion pair. It is related by crossing to the 
gluon helicity-flip skewed pion distribution introduced by Ji and Hoodbhoy 

i- 



Properties of two-pion distribution amplitudes 

Let us now briefly discuss properties of two-pion distribution amplitudes. 
Quark 27rDA's were discussed extensively in Ref. [Q. The gluon distribution 
amplitude $"-^(2;, C, W^) (0) is constrained by soft pion theorem and crossing 
relations. The soft pion theorem for [z^C^W"^) can be obtained in the 
same way as for the quark distribution [Q. It reads: 

$«(^, ( = 0, = 0) = C = 1, _ 0) = . (9) 

Following [Q, we decompose both quark and gluon 27rDA's in conformal 
and partial waves. For quark isosinglet 27rDA the decomposition reads 

00 n+l 

W') = 6z{l - 2) E E BUW)Cl/\2z - l)mC - 1), (10) 

n=l 1=0 
odd even 

while for gluon 27rDA's we have 

00 n+2 

$^(z,c,w^^) = 30 ^2(1-^)2$: Y.^nimci/'{2z-i)mc-i), (n) 

n=0 1=0 
even even 

00 n+2 

<!>^^{z,C,W') = 10z\l-zfY: T.AnFmCl/\2z~l)P/'i2C~l). (12) 

n=0 1=2 
even even 

Here Pi{x) is a Legendre polynomial and Pi"{x) = dP Pi{x) / dx"^ . With this 
choice of the polynomial basis for the (^-expansion the above decomposition 
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is directly related to the partial-wave expansion of the resulting two-pion 
production amplitude. 

Additional constraints for quark and gluon 27rDA are provided by the 
crossing relations between 27rDA's and known parton distributions in the 
pion. The derivation of such relations for quark 27rDA's see ^ . For gluon 
27rDA the derivation is analogous. The crossing relations have the form: 



2 2iV + i /"^ i 

(13) 



^^-2,iv(0) - 5 (jv + i)('^ + 2) Jo '^'^ ''''"'^.{x) , 

where qT^^x), qwix) and g-jri^) are usual quark, antiquark and gluon distribu- 
tions in the pion. Using these relations one can easily derive the normaliza- 
tion of 27rDA's at = 0: 



dz{2z - l)^Q{z, C, W' = 0) = - A M«C(1 - 
JO JMf 

C dz $^(^, C, W^ = 0) = -2 M^C{1 - C) . 
Jo 



(14) 



Here M2 and are momentum fractions carried respectively by quarks 
and gluons in the pion. 

Soft pion theorems and crossing relations constrain ^^{z,(,W^) and 
^'^{z,C,W^) only at = 0. For higher W"^ one can apply dispersion 
relation analysis of Ref. see this paper for more details. 

Under evolution and mix with each other ||^. For asymptotically 
large one obtains (taking into account soft pion theorem (|^)) 

C, w^ = o) = -^^^ z\i - zf C(i - , 

^20 ^ (15) 

C, w' = Q) = ^ ^(1 - {2z - 1) c(i - • 

The distribution amplitude does not mix with quarks under evolution 
and vanishes at asymptotically large . Taking into account contribution 
corresponding to the smallest anomalous dimension leads to the following 
"preasymptotic" expression 



C, W'' = Q)= 30 z\l - zf Al^{Q) , (16) 

where Aq^{0) is unknown constant depending on the normalization point fi. 
This constant is related to matrix element of local operator G^^^ (0)6*^^^(0) 
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in a pion and generically can be of order unity. In eq. ( ]T6| ) dependence on 
the scale is shown exphcitly. The corresponding anomalous dimension 
7-^*^ was calculated in Ref. |^: 



Helicity structure of hard partonic amplitudes 

Let us consider the hard amplitudes 7* + 7 — >■ gg and 7* + 7 — (? + Let 
A*, A denote the respective helicities of the virtual and real photon and Ai, A2 
the helicities of produced partons, quarks or gluons. Due to parity invari- 
ance, the hehcity amplitude MyxMM must obey Ma*a,AiA2 = ^-a*-a,-Ai-A2- 
Because hard scattering occurs coUinearly, conservation of the angular mo- 
mentum along the collision axis results in the helicity conservation A* — A = 
Ai + A2. For the case when photons create a quark- ant iquark pair one finds 
the following amplitudes 

L^ = : M„ 1 _i or M„ _i 1 

2 2'2 ^^^^ 

L, = ±l: Mo_,^i^ior M,,^_.^_i 

Note, however that the = ±1 amplitudes require a helicity flip along the 
fermion line and therefore vanish to the leading twist accuracy. It is easy 
to show that the quark 27rDA (||) describes a subsequent evolution of a 
coUinear quark-antiquark pair with = into two pions. For the case when 
two gluons emerge from the hard collision one finds 

= : Mn,i,-i ot Mu,-i,i 

(19) 

= ±2 : Mi_i,n 

A short calculation shows that the evolution of two collinear gluons with Lz = 
and Lz = ±2 into two pions is described by the distribution amplitudes 
and <l>-^'-^, respectively. Note that the latter results in a nontrivial dependence 
of the cross section on azimuthal angle ip between the lepton scattering plane 
and the pion plane. As it follows, one can extract information about tensor 
gluon distribution ^^'^{z,(,W^) from the corresponding differential cross 
section. 



The leading twist 7*7 irir amplitude to the NLO accuracy 
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Here we present the NLO results for the leading-twist contribution to the 
amplitude @. First, we recall the LO contribution discussed in [|l|. The 
Born graphs contributing to the LO coefficient function are shown in Fig. 3. 



+ crossed 

Fig. 3 LO contributions to the coefficient function 



To this accuracy the colliding photons can produce a quark- ant iquark 
pair only. A direct calculation yields [Q 

T'" = ^ (E ^9 {-g^T [ dz c, W) - (20) 

Note that up to QCD evolution effects the amplitude scales as a function 
of at fixed invariant mass W"^ and C,. One-loop corrections to quark and 
gluon coefficient functions arise from diagrams depicted in Figs. 4 and 5. 




+ crossed 



Fig4. One-loop contributions to the quark coefficient function 



The final result for amplitudes corresponding to the hard subprocess with 
Lz = can be written as: 



Ee9(-OT / dz^'^iz.C.W 



(21) 

As discussed in the previous section, gluons give rise to the hard amplitude 
with Lz = ±2 which results in the contribution 



rpflU 



dz^^''iz,C,W'') ^CUz) 



f22) 



Here K^^K'^ = Kj^K'^ — l/2{—g^'^)±Kj_ and K± is the transverse component 
of the pions momenta (Rf). 
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+ crossed 



Fig. 5 Diagrams for gluon coefficient functions 

Below we present explicit formulae for the coefficient functions Cg,Cg and 
introduced in 



£). Note that they can be obtained either by a direct 
calculation, or in a much simpler way by using the crossing symmetry of 
hard partonic amplitudes which relates them to the previously computed 



NLO coefficient functions for the DVCS process |10, O, O. 



In'fl-^) 



In^z 



+ 



In^;^ 



In^fl 



z(l-z) z{l-z) il-z) 



+3- 



Cl{z) 



Ch 



In z 
1 



Mil-z] 
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+ - - 

z 



l-z) 
9 

(1^ 



^2(1 

+2z{l - 



-zf 

z)\yi[z{1 



z'ln^ z+{l- zf\n^{l- z) + 

z)] - Azlnz - z) ln(l 



(23) 



^)] 



zl 



To get the complete NLO expression one has to include corresponding 
corrections to the and distribution amplitudes. Previously, the NLO 
corrections to the distribution amplitudes in the MS scheme have been found 
numerically small |14|, |T^. In the following numerical analysis we have 
assumed that their effects can be neglected as compared with uncertainties 
introduced by assumptions necessary to model the distribution amplitudes 
themselves. 



Numerical results 

Here we present numerical estimates of the NLO effects in the leading- 
twist amplitude for 7*7 — vrvr process. For this purpose, we assume that 
the corresponding 27rDA's are not far from their asymptotic, respectively 
preasymptotic forms, see eqs. (|T5|,|T6|). To continue 27rDA's from the point 
W"^ = to the higher W"^ values one applies the Watson final state interaction 
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theorem |T6[ and Omnes solution [|T^ to corresponding dispersion relations, 
see Ref. Q for details. With these assumptions, the = amplitude reads: 



-9 



IJ,V\ 



-40 



Nf + ACf 



1 



87 a,(g2 



9 Atx 



X 



/o(W^) Po(cOS^cm) - ^ /2W P2(cOS^e 



• (24) 



Let us note that the asymptotic 27rDA's eqs. (|r3|,|TB|) allow production of 
the pions only in S- and D-waves. Deviations of the 27rDA's from their 
asymptotic form result in appearence of higher waves of produced pions. 

The functions fo,2{W), so-called Omnes functions, can be related to vrvr 
phase shifts 5o(W^) and 52(^) using Watson theorem and dispersion relations 
derived in M: 



fi{W) = exp 



1/1/2 „oo 

i6?(W) + Re / ds-, 

' 71 J Aral s{s-W^-iQ) 



(25) 



IS 



and for the esti- 



The corre- 



The Omnes functions were analyzed in details in Ref. 
mates of the cross sections one can use the results of Ref. 
sponding estimates will be given elsewhere. 

The constant C in eq. (12^) plays a role of an integration constant in the 
Omnes solution to the corresponding dispersion relations. From the soft pion 
theorem it follows that (7 = 1 + O(m^). 

Let us note that the expression ( [25| ) is valid strictly speaking only in the 
elastic region (4m^ < W"^ < 16m^). It should not be complicated to extend 
its range of applicability including the contributions of higher intermediate 
states (probably the most important is the contribution of KK) in the disper- 
sion relations. Note that the accuracy of the Omnes solution can be increased 
if one has an independent information about, say, slopes of functions fi{W) 
at small W. Such information can be provided by effective low energy models 
of QCD |T9[. For example, using the instanton model calculations oi Bni{W) 
at low energies 0] one finds the constant C in eq. (|24D to be equal to: 

C = 1 + 6 +0(m^) with b ^ -1.7 GeV-^ . (26) 

The Omnes functions contain information about the isoscalar vrvr reso- 
nances as well as about the nonresonant background. Due to the resonance 
contributions fo(W) is enhanced around W = 1 GeV (/o(980) resonance) 
whereas f2(W) has a peak at W = 1.275 GeV (/2(1270) resonance). Let us 
note that close to the threshold W = 2m.„ production of pions occurs mostly 
in the 5*— wave. The corresponding Omnes function foiW) is fixed near the 
threshold by the Chiral Perturbation Theory [|18[ to be: 



h{w) = 1 + 



2W^ 



mz 



327r2/2 



P loj 



+ 2 + tTT(3\ + 



1927r2/.2 



(27) 
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Numerically, from equation ( p4D it follows that the NLO correction reduce 
the LO = amplitude by about 30% and that the major NLO correction 
is due to the gluon 27rDA. 

With given by equation ( |I6D the = ±2 amplitude reads 



2Y " 47r 



X 20 f ^4^1 A^m/2W- (2J 



Note that the lowest possible partial wave of produced pions in this case is 
that with 1 = 2. It implies that the Lz = ±2 amplitude is strongly enhanced 
around W = 1.275 GeV due to contribution of /2(1270) resonance. 

Studies of the NLO effects using models of 2'7rDA's beyond the asymp- 
totic approximations ([l5| , p!6D will be published elsewhere (see for model 
predictions for the LO amplitude). 

The interference of = and Lz = ±2 amplitudes produces the az- 
imuthal asymmetries ~ cos 2(p in the cross section of e + 7 — > e' + tttt re- 
actions, where (p is the angle between the lepton scattering plane and the 
pion plane, respectively. As the energy dependence of the vrvr phase shifts 
is known from low-energy experiments, the W^— dependence of various asym- 
metries can be to great extend fixed by the low-energy data. This makes the 
reaction 7*7 vrvr one of the cleanest possible sources of information about 
the tensor gluon distribution amplitude 

The process of the interest 7*7 — >■ tttt in the reaction e -|- 7 — > e' + tttt 
competes with electromagnetic background process when the pion pair is 
created by a bremsstrahlung photon. As the bremsstrahlung results in the 
pion pair being produced in a C— odd state this contribution is absent for 
the TT^TT^ final state. In addition, the interference of the background process 
with 7*7 —>■ 7r+7r~ can be selected by various charge asymmetries 

To our best knowledge at present there are no dedicated measurements 
of the process 7*7 tttt. To the contrary, the processes with real photons 
were studied extensively, see e.g P3 and P3] for tt^tt^ and tt+tt^ channels, 
respectively. Let us note that in the CLEO measurements of 77? transition 
formfactor at large |^ the 7*7 — >■ 7r°7r° were considered background 
events. We strongly recommend to reanalyze this data in order to extract 
the 7*7 — > vt^tt" process. 

To summarize, our estimates suggest that the NLO corrections to the 
leading-twist 7*7 —>■ tttt amplitude are not small and therefore of large im- 
portance for analysis of an experimental data. Moreover, we have shown 
that at the NLO accuracy the amplitude contains a novel tensor gluon dis- 
tribution amplitude, which can be studied via azimuthal asymmetries in the 
cross-section. 
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